Bessel polynomials, double factorials and context-free grammars* 
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\ Abstract 



The purpose of this paper is to show that Bessel polynomials, factorials and Catalan 
triangle can be generated by using context-free grammars. 
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1 Introduction 

The grammatical method was introduced by Chen [2] in the study of exponential structures in 
combinatorics. Let A be an alphabet whose letters are regarded as independent commutative 
iT) ' indeterminates. Following Chen [2], a context-free grammar G over A is defined as a set of 

substitution rules replacing a letter in A by a formal function over A. The formal derivative 

O ! 

£Sj , D is a linear operator defined with respect to a context-free grammar G. For example, if 

G = {x — s> xy, y — > y}, then 

>< 1 D(x) = xy, D(y) = y, D 2 {x) = x{y + y 2 ), D 3 (x) = x(y + 3y 2 + y 3 ). 



For any formal functions u and v, we have 

df(u) 

D(u + v) = D(u)+D(v), D{uv) = D(u)v + uD{v) and D(f(u)) = ^^-D(u), 

ou 

where f{x) is a analytic function. Using Leibniz's formula, we obtain 

D n (uv)=p^D k (u)D n - k (v). (1) 

Let [n] = {1, 2, . . . , n}. The Stirling number of the second kind {^} is the number of ways 
to partition [n] into k blocks. Let S n denote the symmetric group of all permutations of [n]. 
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The Eulerian number (z) enumerates the number of permutations in S n with k descents (i.e., 
i < n, n(i) > ir(i + 1)). The numbers (?) satisfy the recurrence relation 

(iH*+i><r)+<"-Ki"}' 

with initial condition (fy = 1 and boundary conditions (^) = for k > 1. There is a close 
relationship between context-free grammars and combinatorics. The reader is referred to (3J [7] 
for recent results on this subject. Let us now recall two classical results. 

Proposition 1 (J2J Eq. 4.8]). If G = {x — >• xy,y — > y}, then 

= *£{?}**■ 
k=i ^ j 

Proposition 2 (|H Section 2.1]). If G = {x — > xy,y — >• xy}, £/ien 

n-i , > 

fc=0 ^ ' 

The purpose of this paper is to show that Bessel polynomials, factorials and Catalan triangle 
can be generated by using context-free grammars. 



2 Bessel polynomials 

The well known Bessel polynomials y n (x) were introduced by Krall and Frink [5] as the polyno- 
mial solutions of the second-order differential equation 

r,d 2 y n (x) , .dy n (x) , . , . 

x ±y~ + {2x + 2) ^^ = n(n + i)yn{x) - 

The Bessel polynomials y n (x) are a family of orthogonal polynomials and they have been exten- 
sively studied and applied (see [3 A001498]). The polynomials y n {x) can be generated by using 
the Rodrigues formula 

y n (x) = — - — [x ex). 
y K ' 2™ dx n V J 

Explicitly, we have 

k=0 y ' 

These polynomials satisfy the recurrence relation 

y n +i(x) = (2n + l)xy n (x) +y n ^ 1 (x) for n > 0, 

with initial conditions y_i(x) = yo(x) = 1. The first few of the polynomials y n {x) are 

yi{x) = l + x, 

y2(x) = 1 + 3x + 3x 2 , 

y 3 (x) = 1 + 6x + 15x 2 + 15x 3 . 

We present here a grammatical characterization of the Bessel polynomials y n (x). 
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Theorem 3. If G = {a — > ab, b — > b 2 c, c — > be 2 }, then 

D n (ab) = ab n+1 y n (c) for n > 0. 

Proof. Let 

(n + fc)! 

a(n, re) - 



2 fc (n - fc)!fc!" 

Then y n (c) = J2k=o a ( n ' &)c fc - It is easy to verify that 

a(n + 1, re) = a(n, re) + (n + k)a(n, k-1). (2) 

For n > 0, we define 

n 

D n {ab) = ab n+l E(n, k)c k . (3) 

k=0 

Note that D(ab) = ab 2 (l + c). Hence £(1,0) = a(l,0), J B(l, 1) = a(l, 1). It follows from © that 



D n+1 {ab) = D{D n {ab)) = ab n+2 ^ E(n, k)c k + ab n+2 ^(n + re + l)^(n, re)c* 



fc=0 fc=0 

Therefore, 

E(n + 1, fe) = E(n, re) + (n + k)E(n, k - 1). 

Comparing with ([2]), we see that the coefficients E(n, k) satisfy the same recurrence relation and 
initial conditions as a(n,k), so they agree. □ 

For the context-free grammar 

G = {a -> ab, b -> 6 2 c, c -> 6c 2 }, 

in the same way as above we find that 

D n (a 2 b) = 2 n a 2 b n+1 y n (|) for n > 0. 

By Theorem El we obtain D k (a) = ab k yk~\(c) for k > 0. The double factorial of odd numbers 
are defined by 

(2n- 1)!! = 1-3-5 (2n - 1), 

and for even numbers 

(2n)!! = 2-4-6 (2n). 

As usual, set (—1)!! = 0!! = 1. It is clear that 

D n (b) = (2n - l)!!6 n+1 c n for n > 0. 
By ([I]); the following corollary is immediate. 
Corollary 4. For n > 0, we have 



Vn{x) = f>n -2k- 1)!! (?) y k -i(x)x n - 



fe 



3 Polynomials associated with diagonal Pade approximation to 
the exponential function 



The Pade approximations arise naturally in many branches of mathematics and have been ex- 
tensively investigated (see [6] for instance). The diagonal Pade approximation to the exponential 
function e x is the unique rational function 

R (x) - 

Kn[x) - p n (- x y 

where 

P n {x) = Y J M{n,k)x n - k and M(n, fc) = /" + ,^' , . 
£^ [n — ky.kl 

Clearly, P n (l) = y n (2), where y n (x) is the Bessel polynomials. It is easy to verify that the 
numbers M(n, k) satisfy the recurrence relation 

M(n + 1, jfe) = M(n, k) + (2n + 2k)M(n, k - 1). (4) 

The first few of the polynomials P n {x) are given as follows (see [9j Al 13025]): 

P (x) = 1, 

Pi Or) = x + 2, 

P 2 (a;) = x 2 + 6x + 12, 

P 3 (x) = x 3 + 12x 2 + 60x + 120. 

We present here a grammatical characterization of the polynomials P n {x). 
Theorem 5. If G = {a — >■ a& 2 , b — > b 3 c 2 , c — > 6 2 c 3 }, i/ien 

D>6 2 ) = ab 2n+2 c 2n P n (±) ■ 

Proof. For n > 0, we define 

n 

ZJ"(a& 2 ) = a6 2n+2 ^iY(7i,A;)c 2fc . (5) 

fc=0 

Note that D(ab 2 ) = ab 4 (l + 2c 2 ). Hence JV(1,0) = Af(l,0), JV(1, 1) = M(l,l). It follows 
from (J5j) that 



D n+1 (ab 2 ) = ab 2n+4 N(n, k)c 2k + ab 2n+i ^(2n + 2k + 2)iV(n, /fe)c 2 



„2fc+2 

k=0 k=0 

Therefore, 

N(n + 1, jfe) = N(n, k) + (2n + 2k)N(n, k — 1). 

Comparing with ([4]), we see that the coefficients N(n,k) satisfy the same recurrence relation 
and initial conditions as M(n, £;), so they agree. □ 
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Along the same lines, we immediately deduce the following corollary. 
Corollary 6. Let y n {x) be the Bessel polynomials. If G = {a — >• ab 2 , b — » 6 3 c 2 , c — > 6 2 c 3 }, ifoen 



D n (a 2 b 2 ) = 2 n a 2 b 2n+2 y n \ 



c 2 ) 



4 Double factorials 

The following identity was studied systematically by Callan [TJ Section 4.8]: 

J2 m ( 2n ^i x ) ( 2n - 2A: - i) 11 = ( 2n - ( 6 ) 

As pointed out by Callan PQ, the identity ([6]) counts different combinatorial structures, such as 
increasing ordered trees of n edges by outdegree k of the root and the sum of the weights of all 
vertices labeled k at depth n — 1 in the Catalan tree (see [9j A102625]). 
Let 

R{n,k) = kl(^ n ~ k ~ l \{2n-2k-l)\\. 

Thus, 2^fc=i R( n i k) = (2n — 1)!!. It is easy to verify that 

R(n + l,k) = (2n - fc).R(n, fc) + fcl?(n, fc - 1), (7) 

with initial conditions R(0, 0) = 1 and R(0,k) = for fe > 1 or k < 0. For n > 1, let 
R n {x) = Y^!k=i R{ n : k)x k . The first few of the polynomials R n (x) are 

i? 2 (z) = x + 2x 2 , 

i? 3 (a;) = 3x + 6x 2 +6x 3 , 

fl4(x) = 15x + 30x 2 + 36x 3 + 24x 4 . 

Theorem 7. // G = {a — >■ a 2 6, 6 — >• 6 2 c, c — >• 6c 2 }, i/ien 

n 

D n (o) = a6 n R ( n > k)a k c n ~ k for n > 0. (8) 
fc=i 

Proo/. Note that D(a) = o?b and D 2 (a) = ab 2 (ac + 2a 2 ). For n > 1, we define 

n 

D n (a) = ab n r(n, k)a k c n ~ k . (9) 

k=l 

Hence r(l, 1) = R(l, 1), r(2, 1) = R(2, 1) and r(2, 2) = R(2, 2). It follows from Q that 



D(£P(a)) = ao n+1 ^(2n - fc)r(n, A:)a fc c n - fc+1 + ab n+1 J^(k + l)r(n, k)a k+1 c n 



k=l k=0 



Therefore, 

r(n + 1, k) = (2n — k)r(n, k) + kr(n, k — 1). 

Comparing with d7J), we see that the coefficients r(n, k) satisfy the same recurrence relation and 
initial conditions as R(n, k), so they agree. □ 

In the following discussion, we also consider the context-free grammar 

G = {a ->• a 2 b, b ->• 6 2 c, c 6c 2 }. 

Note that 

D(a6) = a 2 b 2 + a6 2 c, L> 2 (a6) = a6 3 (3c 2 + 3ac + 2a 2 ). 

For n > 0, we define 

n 

D n (a6) = a6 n+1 F(n, fc)a fc c n " fc . (10) 

fc=0 

It follows that 



D(D n (ab)) = ab n+2 ^(2n — k + l)H(n, k)a k c n - k+1 + ab n+2 ^(k + l)F(n, k)a k+1 c r ' 



k=0 k=0 

Hence, the numbers H (n, k) satisfy the recurrence relation 



H(n + 1, k) = (2n - k + l)#(rt, it) + kH(n, k - 1), (11) 

with initial conditions #(1,0) = H(l, 1) = 1 and H(l,k) = for fc > 2 or fc < 0. Using ([TT]) . it 
is easy to verify that 

(2n-fc)! 

H(n.k) = r~, ttt- 

v ' > 2 n ~ k (n-k)\ 

It should be noted that the numbers H(n, k) are entries in a double factorial triangle (see [9j 
A193229]). In particular, we have H(n,0) = (2n-l)U,H(n,n) = n\ and Y2=o H ( n i k ) = ( 2n ) !! - 
Moreover, combining dH), ([8]) and (|10p . we obtain 



iT(ra, fc) = (2n - 2m - l)\\R(m, k) 



m=k 

for n > 1 and < A; < n. 
For n > 1, we define 

x(x + 2)(x + 4) • • • (x + 2n - 2) = S^p(n, k)x k 

and 



fe=i 



(a; + l)(x + 3) • • • (x + 2n - 1) = ^ g(n, /c)x fc . 

fc=0 

The the triangular arrays {p(n, fc)}n>i,l<fc<n and {q(n, fc)}n>i,o<fc<n & re both double Pochham- 
mer triangles (see [9l A039683.A028338]). The following theorem is in a sense "dual" to Theo- 
rem and we omit the proof for brevity. 
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Theorem 8. If G = {a — > ab 2 , b — > b 2 c, c — > 6c 2 }, i/ien we /iaue 

n 

D n (a) = ab n J2p( n ^)b k c n - k 
k=l 

and 

n 

D n (ab) = ab n+1 ?(", k)b k c n ~ k . 
k=0 

Set p(0, 0) = ^(0,0) = 1. By ([I]), we immediately obtain 

q(n, k)= ( n )(2n-2m- l)\\p(m, k) 

m=k 

for n > and < k < n. 

5 Catalan triangle 

The classical Catalan triangle is defined by the recurrence relation 

T(n, k) = T(n - 1, k) + T(n, k — 1), 

with initial conditions T(0, 0) = 1 and T(0, fc) = for k > or k < (see A009766]). The 
numbers T(n, fe) are often called ballot numbers. Explicitly, 

m / fn + k\n — k + 1 „ , /N 

T(n, fc) = 7 — for < A; < n. (12) 

\ k J n + 1 

Moreover, Xlfc=o -^( n ' ^) = -^( n + l 5 n + 1) = C( n + 1), where C(n) is the well known Catalan 
number. Catalan numbers appear in a wide range of problems (see [8] for instance). 
It follows from (JT2J) that 

(n + 2)T(n + 1, k) = (n - Jfe + 2)T(n, fc) + (2n + 2k)T(n, k - 1). (13) 

This recurrence relation gives rise to the following result. 

Theorem 9. If G = {a -> a 2 6 2 , 6 ->• 6 3 c 2 , c b 2 c 3 }, then we have 

n 

D n {a 2 b 2 ) = {n + l)\a 2 b 2n+2 T(n, k)a n ~ k c 2k . 

k=0 

Proof. It is easy to verify that D(a 2 b 2 ) = 2a 2 b 4 (a + (?) and D 2 (a 2 b 2 ) = 3!a 2 6 6 (a 2 + 2ac 2 + 2c 4 ). 
For n > 0, we define 



D n (a 2 b 2 ) = (n + l)!a 2 6 2n+2 ^ t(n, fc)a n " fc c 2 



^n—k 2k 

k=0 

Note that 



1 + 1) , Q 2^n+4 = D n " fc + 2 ^ n ' k)a n ~ k+1 c 2k + ]T(2n + 2k + 2)i(n, £>™- fc c 2fc+2 . 



(n + l)._ 

v ; fc=0 fc=0 
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Thus, we get 



(n + 2)t(n + l,k) = (n-k + 2)t(n, k) + (2n + 2k)t(n, k - 1). 

Comparing with (|13p . we see that the coefficients t(n,k) satisfy the same recurrence relation 
and initial conditions as T(n, k), so they agree. □ 

In the same way as above we find that if G = {a — > a 2 b 2 , b — > b 3 c 2 , c — > 6 2 c 3 }, then 

n 

D n {ab 2 ) = n\ab 2n+2 Y J 

k=0 

and 

n-l 

D n {b) = Y[{Ak + l)b 2n+1 c 2n . 

k=0 

It should be noted that ( n ^ fc ) is the number of lattice paths from (0, 0) to (n, k) using steps 
(1,0) and (0,1) (see A046899]) and rifc=o( 4fc + x ) is the Quartic factorial number (see 
A007696]). 
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